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Abstract 

We present a measurement of the unitarity triangle angle 03 using a Dalitz plot analysis of 
the three-body decay of the neutral D meson from the — ► DK*^ process. Using a 253 fb^ 1 
data sample collected by the Belle experiment, we obtain 56 signal candidates for B^ — > DK*^ 
where the neutral D meson decays into Kstt + tt~. From a maximum likelihood fit we obtain 
03 = 112° ± 35°(stat) ± 9°(syst) ± 11° (model) ± 8° (nonresonant B ± -» DKs^). 

PACS numbers: 13.25.Hw, 14.40.Nd 
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INTRODUCTION 



Determinations of the Cabbibo-Kobayashi-Maskawa (CKM) [Jj matrix elements provide 
important checks on the consistency of the Standard Model and ways to search for new 
ohysics. Various methods using CP violation in B — > DK decays have been proposed 
l H El El to measure the unitarity triangle angle 4> 3 . The sensitivity to the angle </> 3 comes 
:'rom the interference of two amplitudes producing opposite flavors of neutral D meson. Such 
a mixed state of neutral D will be called D. For example, in the case of B + decay, the mixed 
state is D + = D° + re ie +D°. Here r is the ratio of the suppressed and favored amplitudes (it 
is expected to be of order of 0.1-0.3), the total phase between the D° and D° is 9 + = <p 3 + 5, 
where 5 is the strong phase difference between suppressed and favored decays. Analogously, 
for the decay of B~, one can write £)_ = D° + re ie ~ D° with = —6 3 + 5. 

Three body final states common to D° and D°, such as Ks7r + 7V~ j6|, have been suggested 
as promising modes for the extraction of 3 . The Dalitz plot density of D gives immediate 
information about r and 9±, once the amplitude of the D° decay is known. The amplitude 
of the D + decay as a function of Dalitz plot variables m 2 + = m 2 KsiT+ and m 2 _ = m 2 KsiT _ is 

M + = f(m 2 + ,m 2 _) + re i,h+iS f(m 2 _, m 2 + ), 
where f(m+, m 2 _) is an amplitude of the D° —>■ As7r + 7r~ decay. Similarly, 

M_ = f(m 2 _,m 2 + ) + re-^+^fiml, m 2 



The D° — > Ks7r + 7r - decay model can be determined from a large sample of flavor-tagged 
D° — > A'57r + 7r~ decays produced in continuum e + e~ annihilation. Once that is known, a 
simultaneous fit of B + and B~ data allows to separate the contributions of r, 03 and 8. 
Refer to 0, 0] for a more detailed description of the technique. 

The method described can be applied to other modes as well as B^ — > DK ± decay. 
Excited states of neutral D and K mesons can also be used, although the values of 6 
and r can differ for these decays. Previously Belle 0, and Babar collaborations have 
performed the analyses using this technique involving B^ 1 — > DK ± and B^ — > D*K ± decays. 
This paper describes the analysis using B ± — > DK*^ 1 decay. 



EVENT SELECTION 

We use a 253 fb _1 data sample, corresponding to 275 x 10 6 BB pairs, collected by the 
Belle detector 0. The decay chain B ± -> DK** with D -> K s it + ti~ and K*± -> A s vr ± is 
selected for the analysis. 

The requirements on the quality of the charged tracks are the same as in our previous 
analysis |8j. To select A* ± mesons, we require the invariant mass of Kgir^ to be within 
50 MeV/c 2 of its nominal mass. For the selection of neutral D meson, the invariant mass 
Mk s tttt of its decay products is required to satisfy \M Ks7rn — m D \ < 15 MeV/c 2 . 

The selection of B candidates is based on the center-of-mass (CM) energy difference 
AE = — E hcam and the beam-constrained B meson mass M bc = \/E 2 c&m — {^Pi) 2 , 

where i?beam is the CM beam energy, and Ei and Pi are the CM energies and momenta of the 
B candidate decay products. The requirements for signal candidates are 5.272 GeV/c 2 < 
M bc < 5.288 GeV/c 2 and \AE\ < 0.022 GeV. 
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To suppress background from e + e~ — > qq (q = u, d, s, c) continuum events, we require 
| cos^thrl < 0.8, where #thr is the angle between the thrust axis of the B candidate daughters 
and that of the rest of the event. We also require the helicity angle of K*^ to be | cos#hci| > 
0.4. For additional background rejection, we use a Fisher discriminant composed of 11 
parameters [ill ]: the production angle of the B candidate, the angle of the B thrust axis 
relative to the beam axis and nine parameters representing the momentum flow in the event 
relative to the B thrust axis in the CM frame. We apply a requirement on the Fisher 
discriminant that retains 95% of the signal and rejects 30% of the remaining continuum 
background. 

The B ± -> DK*± selection efficiency (4.1%) is determined from a Monte Carlo (MC) 
simulation. The number of events passing all selection criteria is 56. The AE and M^ c 
distributions for B ± — > DK*^ candidates are shown in Fig. ^ The background fraction 
is determined from a binned fit to the AE distribution, in which the signal is represented 
by a Gaussian distribution with mean AE = and the background is modeled by a linear 
function. The contributions in the signal region are found to be 36 ± 7 signal events and 
13 ± 2 background events. The overall background fraction is 27 ± 5%. 
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FIG. 1: (a) AE and (b) Mb c distributions for the i? — > DK*^ candidates. Vertical lines show 
the signal region. The histogram shows the data; the smooth curves are the fit result. 



DALITZ PLOT ANALYSIS OF B ± -» DK*± DECAY 

The Dalitz plot distributions for the D — > Ks^ + ^ decay from B^ — > DK*^ are shown 
in Fig. El 

The Dalitz plot fit procedure is similar to that used in our previous analysis 0. The 
model of D° —* Ks^ + 7r~ decay is based on the 253 fb -1 sample of continuum D*^ — > Dn^ 
decays and is also the same as in our B± -> D(*)K ± analysis. 
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FIG. 2: Dalitz plots of D -> KsTr+vr" decay from (a) 5+ -> and (b) 5~ -> DK*~. 



TABLE I: Fractions of different background sources. 



Background source 


Fraction 


qq combinatorial 

BB events other than B ± -> DK*^ 
Combinatorics in D° decay 


18.0 ± 4.5% 
9.0 ± 1.5% 
0.6 ± 0.1% 


Total 


27 ± 5% 



We consider three sources of background (see Table U}, an d determine the fraction and 
Dalitz plot shape for each component. The largest contribution comes from two kinds 
of continuum events: random combinations of tracks, and correctly reconstructed neutral 
D mesons combined with a random K$ and a pion (or a random K*^). The fraction of 
candidates for these events is estimated to be 18.0 ± 4.5% using an event sample in which 
we make requirements that primarily select continuum events but reject BB events. The 
shapes of their Dalitz plot distributions are parameterized by a third-order polynomial in 
the variables m 2 + and m 2 _ for the combinatorial background component and a sum of D° 
and D° shapes for real neutral D mesons combined with random K$ and a pion. 

The background from BB events is subdivided into two categories. The DK*± combi- 
nations coming from the decay of D^*> from one B meson and K$ and ^ from the other 
B decay constitute the largest part of the BB background. The fraction of this source is 
estimated to be 9.0 ± 1.5% using a MC study. The B ± — > DK*^ events where one of the 
neutral D meson decay products is combined with a random kaon or pion were studied 
using a signal MC data set. The estimated background fraction is 0.6 ± 0.1%. The Dalitz 
plot shape for both of these backgrounds is parameterized by a sum of linear functions of 
variables m 2 + and m 2 _ and a D° amplitude. 
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RESULTS 



Fig. 13 shows the constraints on the complex amplitude ratio re lB separately for B + and 
B~ samples. The global minima of the likelihood function are denoted by the crosses. 
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FIG. 3: Constraint plot of the complex amplitude ratio re td for B^ — > DK*^ decay. Contours 
indicate integer multiples of the standard deviation. Dotted contours are from B + data, dashed 
contours are from B~ data. 

A combined unbinned maximum likelihood fit to the B + and B~ samples with r, 3 
and 5 as free parameters yields the following values: r = 0.37 ± 0.18, 03 = 112° ± 32°, 
6 = 353° ±32°. The errors quoted here are obtained from the likelihood fit. These errors are 
a good representation of the statistical uncertainties for a Gaussian likelihood distribution, 
however in our case the distributions are highly non-Gaussian. In addition, the errors for 
the strong and weak phases depend on the values of the amplitude ratio r (e.g. for r = 
there is no sensitivity to the phases). 

As in our 

£± £>(*)x± analysis, we use a frequentist technique to evaluate the statistical 
significance of the 03 measurement and to correct for the bias of the fit procedure. We use 
toy MC pseudo-experiments to obtain the probability density function (PDF) of the fitted 
parameters as a function of the true parameters, followed by a confidence level calculation 
using the Neyman method. The confidence regions for the pairs of parameters (03, t) and 
(03, 5) are shown in Fig.0] They are the projections of the corresponding confidence regions 
in the three-dimensional parameter space. We show the 20%, 74% and 97% confidence level 
regions, which correspond to one, two, and three standard deviations for a three-dimensional 
Gaussian distribution. While the 03 and S values that are determined from the toy MC are 
consistent with those that are determined in the unbinned maximum likelihood fits for both 
DK ± and D*K ± , the corresponding r values are significantly different. This is caused by 
a bias in the unbinned maximum likelihood. Since r is a positive-definite quantity, the fit 
tends to return a larger value for r than its true value, particularly when r is small. 

There are potential sources of systematic error such as uncertainties in the background 
Dalitz plot density, efficiency variations over the phase space, m^ n resolution, and possible 
fit biases. These are listed in Table |TI] The effect of background Dalitz plot density is 
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FIG. 4: Confidence regions for the pairs of parameters (a) (r, ^3) and (b) (03,5) for the 
DK*^ 1 sample. 



TABLE II: Contributions to the experimental systematic error. 



Source 


Ar 


A0 3 (°) 


A<5(°) 


Background shape 


0.027 


5.7 


4.1 


Background fraction 


0.006 


0.2 


1.0 


Efficiency shape 


0.012 


4.9 


2.4 


resolution 


0.002 


0.3 


0.3 


Control sample bias 


0.004 


10.2 


10.2 


Total 


0.03 


13 


11 



estimated by extracting the background shape parameters from the sidebands and by 
using a flat background distribution. The maximum deviation of the fit parameters from 
the "standard" background parameterization is assigned as the corresponding systematic 
error. The effect of the uncertainty in the background fraction is studied by varying the 
background fraction by one standard deviation. The efficiency shape and resolution are 
extracted from the MC simulations. To estimate their contributions to the systematic error, 
we repeat the fit using a flat efficiency and a fit model that does not take the resolution 
into account, respectively. There is no obvious control sample for the mode — > DK*^ 
(with K*± -> Ks^) like B ± -> D^ir* in the case of B ± -> D^K* modes, but we include 
the effect of the possible bias of the B^ — > D-k^ control sample into the systematic error to 
account for possible deficiency in the D° model, or other fit biases. 

The uncertainty due to the model of D° decay is taken to be the same as for the B — > 
processes. However, the analysis of B ± DK*^ has an additional uncertainty due 
to the possible presence of the nonresonant B^ — > DKsff component, which can also be 
considered as a model uncertainty. Since the nonresonant decay is described by the same 
set of diagrams as B ± — > DK*^ a similar CP violating effect should take place but, in 
general, with different r and 5 from the resonant mode. Thus, for the <f>3 measurement 
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using — > DK*^ mode alone, the contribution of — > DKgTT decay can bias the fit 
parameters. To estimate the corresponding systematic uncertainty, we set an upper limit on 
the B ± — > DKsir fraction, and perform a toy MC simulation with nonresonant contribution 
added to determine the fit bias. 

To limit the fraction of the nonresonant component we search for a B^ — ► DKg^ signal 
in the sidebands of the K* invariant mass distribution. Since a specific final state of the 
neutral D meson is not needed for this study, we combine several D decay modes (K7T, 
KsTrTr, Kirn , Ktytity) to increase the statistics. No requirements on event shape parameters 
or K* helicity are applied. This allows to obtain a tighter upper limit on the £? — > DKg^ 
contribution due to higher signal efficiency. 

The following ranges of K* invariant mass m^ s7r are analysed: |m^ s7r — Mk*\ < 50 
MeV/c 2 for B ± -> DK** events and 80 MeV/c 2 < \m Ks7T - M K *\ < 300 MeV/c 2 for 
B^ — > DRgn . The number of signal events is extracted from the fit to AE distribution. 
The fit function is a combination of a Gaussian peak for the signal and a linear function for 
the background. The number of events in the K* sideband is 3 ± 40. The number of events 
in the K* signal region is 321 ±31. Taking into account the shape of Ksir invariant mass 
distribution, we obtain a B^ — » DKg'K fraction in our signal region of 0.3 ± 3.7%. The 
95% confidence level (CL) upper limit on this fraction is 6.3%, which we use for estimation 
of the systematic error due to the nonresonant contribution. 

To study this effect, the Dalitz plot of the D decay is generated according to a decay 
amplitude with a component having different values of r and 5, corresponding to the B^ — > 
DKg7T contribution, added coherently. The resulting Dalitz plot distribution is then fitted 
with the standard technique, and the biases of the fit parameters r, 3 and 5 are obtained. 
The fits are performed for different values of r and 8 for the nonresonant component and for 
different values of the relative phase between B^ — > DK*^ and B ± — > DKs^ amplitudes. 
The maximum bias of the fit parameters is taken as a corresponding systematic error. We 
obtain the following estimates of the uncertainty due to the B — > DKg^ contribution: 
Ar = 0.084, A03 = 8.3°, AS = 49.3°. We quote these errors separately from the other 
systematic uncertainties. The 03 bias is significantly smaller than that for the strong phase 
5, since 03 is obtained from a difference of the total phases for B + and B~ decays, and a 
part of the bias cancels in this case. 

For the final results, we use the central values that are obtained by maximizing the PDF 
of the fitted parameters and the statistical errors corresponding to the 20% confidence region 
(one standard deviation). Of the two possible solutions (03, 5 and 03 + 180°, 5 + 180°) we 
choose the one with 0° < 3 < 180°. The final results are 

r = 0.25^;^±0.09±0.04±0.08, 3 = 112 ±35 ±9 ±11°±8 , 5 = 353°±35 ±8 ±21 ±49°. 

The first error is statistical, the second is experimental systematics, the third is model 
uncertainty and the fourth is the error due to possible S ± — > DKs^ contribution. The 
two standard deviation interval including the systematic and model uncertainties is 34° < 
03 < 190°. The statistical significance of CP violation is 63%. 

CONCLUSION 

We report results of a measurement of the unitarity triangle angle 03 that uses a method 
based on a Dalitz plot analysis of the three-body D decay in the process B ± — > DK*^. The 
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measurement is based on a 253 fb 1 data sample collected by the Belle detector. We obtain 
the value of 3 = 112° ± 35° ± 9° ± 11° ± 8° (solution with 0° < 3 < 180°). The first error 
is statistical, the second is experimental systematics, the third is model uncertainty and the 
fourth is the error due to possible — > DRg-K contribution. The statistical significance 
of CP violation for the combined measurement is 63%. The method allows us to obtain a 
value of the amplitude ratio r, which can be used in other 3 measurements. We obtain 
r = 0.25±°;^ ± 0.09 ± 0.04 ± 0.08. 
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